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ABSTRACT 


Determining the best base stationing for military units can be modeled as a 
capacitated facility location problem with sole sourcing and multiple resource 
categories. Computational experience suggests that cross decomposition, a 
unification of Benders Decomposition and Lagrangean relaxation, is superior to other — 
contemporary methods for solving capacitated facility location problems. Recent 
research extends cross decomposition to pure integer programming problems with 
explicit application to capacitated facility location problems with sole sourcing; 
however, this research offers no computational experience. This thesis implements 
two cross decomposition algorithms for the capacitated facility location problem with 
sole sourcing and compares these decomposition algorithms with branch and bound 
methods. For some problems tested, cross decomposition obtains better solutions in 
less time; however, cross decomposition does not always perform better than branch 
and bound due to the time required to obtain the cross decomposition bound that is 


theoretically superior to other decompositon bounds. 
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EXECUTIVE SUMMARY 

The problem of determining the best stationing plan for military units can be 
modeled as a capacitated facility location problem with sole sourcing and multiple 
resource categories. The capacitated facility location problem accomplishes two 
goals. It determines which facilities to open and how to allocate customer demand 
among the open facilities. For unit stationing this equates to which bases to retain and 
which units to assign to bases. Sole sourcing requires that all customer demand be 
allocated to one facility. For the base stationing application this requirement means 
that a unit is stationed at a single base. In capacitated facility location problems with 
multiple resource categories the customers require more than one type of resource 
from the facilities. For unit stationing problems the resources include those that are 
consumed by the presence of the unit, such as barracks and motor pools, and those 
that are shared by units, such as ranges. 

The most common techniques for solving models like the capacitated facility 
location problem with sole sourcing are branch and bound methods. Some problem 
instances of the capacitated facility location model with sole sourcing require many 
hours to solve using branch and bound. For some instances it is not practical to solve 
the problem optimally, but it is possible to find a feasible solution and a lower bound 
on the cost of the optimal solution. Decomposition algorithms are another common 
technique for solving capacitated facility location problems with sole sourcing. These 
techniques decompose a large or difficult problem into a series of easier problems. 
Recent research has unified the two predominant decomposition approaches, Benders 
decomposition and Lagrangean relaxation, into a decomposition framework called 
cross decomposition. Computational experience has shown that cross decomposition 
can efficiently solve the capacitated facility location problem. This research does not, 
however, provide computational experience with the capacitated facility location 


problem with sole sourcing. 
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This thesis implements two cross decomposition algorithms for the capacitated 
facility location problem with sole sourcing and compares these decomposition 
algorithms with branch and bound methods. These algorithms are implemented using 
a commercial solver on a high speed digital computer. For some problems tested, 
cross decomposition obtains better solutions in less time; however, cross 
decomposition does not always perform better than branch and bound due to the time 
required to obtain the cross decomposition bound that is theoretically superior to 
other decompositon bounds. The favorable results obtained using cross 


decomposition indicate that the method is worthy of further research. 
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I. INTRODUCTION 

This thesis implements two cross decomposition algorithms [Holmberg 1994] for a 
capacitated facility location problem with sole sourcing, compares these solution methods 
with direct solution using branch and bound, and applies a decomposition algorithm to a unit 
stationing model [Dell et al. 1994]. This research evaluates the ability of a cross 
decomposition algorithm to determine lower bounds for the optimal objective function value 
of the capacitated facility location problem with sole sourcing and recommends a method to 
obtain feasible solutions from the cross decomposition algorithm. 

The remainder of this chapter introduces the stationing model, the cross 
decomposition algorithm and the research approach. Chapter Il describes the capacitated 
facility location problem and the cross decomposition method. Chapter III describes the 
application of two cross decomposition algorithms to the capacitated facility location problem 
with sole sourcing. Chapter IV presents and analyzes the test results. Chapter V offers 
conclusions and recommendations for further study. 

A. BACKGROUND 

1. Optimal Stationing of Units at Bases (OSUB) 

The Defense Base Closure and Realignment Act of 1990 [BRAC 1990] 
specifies a process by which major Department of Defense installations can be closed or 
realigned. This law and subsequent legislative amendments require that installations be 
evaluated primarily by military value and cost. To assist the Army with stationing decisions, 
the OSUB (Optimal Stationing of Units to Bases) model [Dell e¢. al. 1994] was developed 
at the Naval Postgraduate School. This model is an elastic mixed integer programming model 
with two criteria: military value and cost. OSUB has recommended closure and realignment 
for Army maneuver and training bases and stationing Army units. OSUB is currently 
implemented with the General Algebraic Modeling System, GAMS, [Brooke et. al. 1992] and 
solved using the integer linear program solver XA [Sunset Technology 1987]. OSUB is a 
capacitated facility location problem with sole sourcing and some special constraints. The 


OSUB model formulation is presented below. 








Indices: 


Data: 





i,i’ bases; 
j units; 


k resources (includes total maneuver acres, contiguous maneuver acres 
(ca), housing, facilities and ranges). 


S, set of units that are currently stationed at base /; 
futil, fixed utility of base 7, 
feost, fixed cost associated with keeping base 7 open; 


vutil,, the difference in variable utility when unit 7 moves to base i (a 
positive difference is a desirable change); 


vcost,, the difference in variable cost when unit 7 moves to base / (a 
positive difference is a higher cost); 


pen,, penalty per unit deviation from resource k at base i (any deviation not 
associated with a military value objective has a penalty of zero); 


co,, operating cost associated with deviating from resource k at base / (any 
deviation not associated with the cost objective has a cost of zero); 


cc, construction cost associated with deviating from resource k at base i 
(any deviation not associated with a construction cost has a cost of zero); 


cap,, capacity of resource k at base j (current stationed unit use is 
subtracted from the capacity for all k # ca); 


r,, resource & utilization by unit j; 
cm,, cost to move unit j to base /; 


maxm the maximum movement cost, 


maxc the maximum one-time realignment cost. 








Binary Variables: 


* close,= 1 if base i is closed and 0 if it remains open, 


* move,,= 1 if unit j moves to base 7. (This variable is defined only for units 


j not already assigned to base .) 


Continuous Variables: 


* e., deviation from resource capacity k at base 7. 


Formulation: 


(1) maximize Z, = pa futil; -d— close; +>, > vutil,, “move; 2d, pen; i, 


t jes, 


(2) minimize Z, “- fcost, -(1—close; +>, > veost; ; ‘move; ; + dd a ene 


(3) 


(4) 
(5) 


(6) 


(7) 
(8) 


(9) 


(10) move, 


i jes, 

subject to: 
> move, .<1 WijeS, 

ty t 
i'#i 
move; , <(1—close;) Vi,j €5; 
move, .2close, Vi,j €S; 

tyJj 
i'#i 
Dig “Move; ; Os, rj p Move, , SCOD,, + ey 
JES; jes, ii 
r;, Move; , <cap;,+€, Vij €S,k=ca 


S| > cm, ;-move; ; < maxm 
ioJ 

Sd icein ein S marc 
ik 


j ={0, 1} Vi,j éS, close, {0,1} Vi 


These objectives and constraints respectively: 


Vi,k #Ca 


€;, 20 Vi,k 


(1) express a comparative measure of military value for units assigned to bases, 


(2) express the cost of unit stationing, 














(3) ensure a unit moves at most once or not at all, 

(4) & (5) ensure that a base closes only if no new units are stationed there and all old 
units move to another base, 

(6) capacitate housing, total maneuver acres, facilities and ranges or measure excess 
demand for such, 

(7) capacitate contiguous maneuver acres or measure excess demand for such, 

(8) limit the maximum unit movement cost, and 

(9) limit the one-time construction cost for realignment. 

Zi Cross Decomposition 

Van Roy [1983, 1986] develops a cross decomposition algorithm that unifies 
Benders decomposition and Lagrangean relaxation and applies it to the capacitated facility 
location problem. This method simultaneously isolates and exploits the primal and the dual 
structure of the problem by successively solving a transportation problem and a simple plant 
location problem. Van Roy presents evidence that this algorithm is superior to several other 
methods and shows that an implementation of the algorithm solves sample problems about 
ten times faster than other methods available at that time. Van Roy also claims that in all 
problems tested his method obtains in just a few decomposition iterations tight lower and 
upper bounds that differ by no more than 0.5%. 

Holmberg [1990, 1994] generalizes the concept of cross decomposition to 
pure integer problems and studies the lower bounds on the optimal objective function value 
of pure integer programming problems. He proves that generalized Benders decomposition 
and generalized cross decomposition yield the best lower bound for pure integer problems. 
He also shows that while either ordinary Benders decomposition or Lagrangean relaxation 
may yield the best lower bound for a particular problem, cross decomposition can 
automatically yield the best of these bounds. He formulates the decomposition problems 
required to apply generalized cross decomposition to a capacitated facility location problem 
with sole sourcing. He provides no computational evidence of the efficiency of these 


methods. 





B. OVERVIEW 
This thesis presents three accomplishments. 
I. Implementing Cross Decomposition 
Two cross decomposition algorithms are implemented for the capacitated 
facility location problem with sole sourcing. The algorithms are coded in C on an IBM 
System 6000 Model 590H using the CPLEX Callable Library. CPLEX [1994] is an 
optimization tool for solving linear and mixed integer optimization problems. The CPLEX 
Callable Library is an object-oriented C library. CPLEX allows the user to build applications 
which solve, modify, and interpret the results of linear and mixed integer programs. To code 
the algorithms we develop techniques to successively generate and solve the decomposition 
problems. 
2; Analysis of Results 
The cross decomposition algorithms are tested and analyzed. We perform 
computational tests using the decomposition algorithms and the CPLEX Mixed Integer Solver 
on problems of various sizes. We compare the results and characterize the performance of 
the decomposition algorithms. This analysis evaluates the capability of the cross 
decomposition algorithms to determine lower bounds on the optimal objective function value 
for capacitated facility location problems with sole sourcing. We also analyze the suitability 
of the algorithms to obtain feasible solutions using the facility configurations and lower 
bounds from the decomposition algorithms. 
3. Application to Unit Stationing Models 
A decomposition approach is selected and used on a unit stationing model 
similar to the OSUB model. The unit stationing model is simple extension of the capacitated 


facility location problem with sole sourcing. 








Il. CAPACITATED FACILITY LOCATION AND CROSS DECOMPOSITION 


A. CAPACITATED FACILITY LOCATION PROBLEMS 

The CFLP describes a wide variety of planning problems. Applications beyond 
capacitated facility location include: lot sizing decisions in production planning; 
telecommunications network design; machine replacement; vehicle routing when capacities 
are not equal [Cornuejols ef. al. 1991]; the stochastic transportation problem; and discrete 


network design [Holmberg 1990]. The mixed integer formulation of the CFLP is presented 


below. 
Indices: 
* 7 facilities; 
* j demand points. 
Data: 
* f, fixed cost to operate facility 7; 
* c, , cost to supply all demand at / by facility 7; 
* d,total demand at /. 
Binary Variables: 


* y= 1 if facility 7 is open and 0 if it is closed. 


Continuous Variables: 


* x, , the proportion of the total demand at j (d) provided by facility 7. 











Formulation: 


minimize SA it Cay Xj 
i ij 
yx 


(CFLP) 


subject to: 


(CFLP 1) ye =] [Aj] Vj 
j 


(CFLP 3) xj SY; [v; ;] Vi,J 


(CFLP 4) x; 2 QO Vi,J y; €{0,1} Vi 


The objective is to minimize total cost, and contains two distinct sets of decision 
variables. The first set of binary decision variables (,) determine which facilities to open. 
The second set of continuous decision variables (x, ,) allocate customer demand to the open 
facilities. The variables in brackets are the corresponding dual variables for the given set of 
constraints. Constraints (CFLP 1) ensure that all demand is met. Constraints (CFLP 2) 
enforce the capacity limits for the facilities. Constraints (CFLP 3) are variable upper bounds 
on the allocation of demand. This third set of constraints is redundant in binary variables, but 
yields a much tighter linear programming relaxation [Van Roy 1986]. 

B. SOLVING THE CFLP 
L; Benders Decomposition 
- Benders decomposition (e.g., [Nemhauser & Wolsey 1988] and [Magnanti & 
Wong 1990]) is a primal solution method. It isolates special structure in the problem by 
fixing primal variables. Benders decomposition is an exact method that solves the CFLP 
optimally by iterating between the primal subproblem and the primal master problem 
described below. Appendix A contains a derivation of the master problem and subproblem. 


The description below provides implementation details. 








a. The Primal Subproblem (PS) 


minimize > Si ‘y+ , on Xi ij 
(PS) i i,j 
j 


subject to: 


(PS 1) > *1j =! [A] W 


i 


(PS 2) Sod, Xi ij <5; Yj Lz; ] Vi 
j 


(PS 3) x; SY; [v; ;] Vi, J 


(PS 4) Xij >0 WVi,j 


The primal subproblem is the linear program obtained when the facility 
configuration is fixed at in the CFLP where y, indicates that y,is fixed at either one or zero. 


PS is a restriction of the CFLP that consequently provides an upper bound on its optimal 


objective function value. The objective function value of the dual of PS is 
» 4; =e Yi -> 3 Vii 
j i i,j 
To maintain consistency with the objective function of the CFLP, we add the same constant, 


(the fixed facility cost), that is added to PS and obtain 


bas: Yi +A, a Yi; ->), Vij 
i J ij 


i 


By duality, the value of the above expression for the optimal solution to PS is the maximum 


feasible value for the given facility configuration. This expression provides the basis for a 








primal cut. The optimal objective function value of the primal master problem must be less 


than or equal to 
(f) (t) (t) 
DAP UG HA 8 MY: 
j i j 


for any set of feasible facility configurations (where the superscript (#) on the dual variables 
is the iteration number). 


b. The Primal Master Problem (PM) 


minimize 
(PM) | 
yp 


subject to: 


(PM 1) p> yO 4G w 8-Day y: Vt 
j i 


j 
(PM 2) Da -y; > od, 
j 


i 


(PM 3) y; €{0,1} Vi 


The primal master problem is obtained by adding a primal cut (PM 1) 
each time the primal subproblem is solved. PM fixes the facility location variables for the 
primal subproblem. PM provides a feasible facility configuration because constraint (PM 2) 
ensures total supply meets total demand. 

When using only a subset of all possible cuts (a relaxation), PM 
provides a lower bound on the optimal objective function value of the CFLP. Since PS 


determines the maximum value of 


¥ 4) GaP DY 
j i j 


10 








for the facility configuration provided by PM, the optimal objective function value of the 
CFLP is identified when the objective function value of PS equals the objective function value 
of PM. 
2; Lagrangean Relaxation 

Lagrangean relaxation [e.g., Nemhauser and Wolsey 1988] is a dual solution 
method. It isolates special structure in the problem by moving complicating constraints into 
the objective function and penalizing complicating constraint infeasibility with a Lagrangean 
multiplier. Lagrangean relaxation solves a relaxation of the CFLP by iterating between the 
subproblem and the master problem described below. 


a. The Dual Subproblem (DS) 


minimize . 
@s) DA Yi +o eo A => » 4 de a a 
iti 


x,y i i,j 


subject to: 


s 1) er =1 W 
i 


oe Yao De 
i J 


(DS 3) Xi j <y, Vij: 


@s4 x; ,2 QO Vi,j y, €{0,1} Vi 


The dual subproblem is the Lagrangean relaxation of the CFLP with 
respect to the total supply constraint (CFLP 2) and the addition of a constraint (DS 2) that 
ensures total supply meets total demand. DS is a relaxation of the CFLP that consequently 
provides a lower bound on its optimal objective function value. DS provides for any set of 


_ dual values fixed facility locations (y,) and customer assignments (x, ,). 


I] 








b. The Dual Master Problem (DM) 


maximize ro) 
(DM) 
[1,0 


subject to: 


(t) (t) (t) (t) 
i i,j i 


j 


M2 ft; 20 Vi 


The dual master problem is obtained by adding a dual cut (DM 1) each 
- time the dual subproblem is solved. With all possible cuts, DM is the Lagrangean dual. DM 
fixes the dual variables for the dual subproblem. DM provides an upper bound on the optimal 
objective function value of the Lagrangean relaxation of the CFLP because it is a relaxation 
of the Lagrangean dual. 
3. Cross Decomposition 

Cross decomposition unifies Benders decomposition and Lagrangean 
relaxation. Figure 1 below motivates the cross decomposition method. Cross decomposition 
iterates in the subproblem phase between the restricted primal subproblem and the relaxed 
dual subproblem described above. These problems successively provide an upper and lower 
bound on the optimal objective function value of the CFLP. These problems may provide 
tight bounds in the subproblem phase, but neither convergence, nor monotonic improvement 
is guaranteed. The convergence tests described below are used to determine when the 
subproblems fail to make progress toward an optimal solution. When a convergence test fails, 
cross decomposition solves a master problem that is formed using the cuts generated in the 


subproblem phase and then continues this phase with the next subproblem. 
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Cross Decomposition 


Primal 
Subproblem 


Dual 
Subproblem 


Primal Master Primal 


Dual Master Problem Convergence 
Problem Test 


Primal 
Subproblem 





Figure 1. 

Cross Decomposition iterates between the primal and dual subproblems in the "subproblem 
phase." It solves a primal or dual master problem when a convergence test fails and then restarts 
the subproblem phase. The algorithm terminates when the objective function values of selected 


problems converge to a lower bound on the optimal objective function value. 


a. The Primal Convergence Test 


The primal convergence test uses the following cuts: 


“ 


t) (t) (t)| « 
>A; +), fi - 2; 5;- Vi -y, < Upper Bound 
j i j , 


If these cuts are satisfied for all ¢ then cross decomposition continues by solving the primal 
subproblem. The primal convergence test uses the cuts in the primal master problem to 
determine if the upper bound can be improved. If any cut ¢ 1s not satisfied, a master problem 


is solved. The primal convergence test is a necessary condition for improving the current 
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upper bound on the optimal obj ective function value. The primal convergence test is 

sufficient to show that either: the primal subproblem can improve the upper bound; or the 

primal subproblem can generate a new cut for the primal master problem [Holmberg 1990]. 
b. The Dual Convergence Test 


The dual convergence test uses the following cuts: 


(t) (t) (t) (ft) | « 
D694 QU Fy! ‘ S 4X55 - 85°93 - Lt, = Lower Bound 
i,j i lj 


If these cuts are all satisfied then cross decomposition continues by solving the dual 
subproblem. The dual convergence test uses the cuts in the dual master problem to determine 
if the lower bound can be improved. If any of these cuts are not satisfied, cross 
decomposition solves a master problem. The dual convergence test is a necessary condition 
for improving the current lower bound on the optimal objective function. The dual 
convergence test is sufficient to show that either: the dual subproblem can improve the lower 
bound; or the dual subproblem can generate a new cut for the dual master problem 


{Holmberg 1990]. 


14 


Ot. CAPACITATED FACILITY LOCATION WITH SOLE SOURCING 


The capacitated facility location problem with sole sourcing (CFLPSS) is an important 


variant of the CFLP model. In this problem all customer demand must be assigned to one 





facility and the second set of decision variables (x, ,) are also binary. The CFLPSS is a pure 
integer linear program with significantly more binary variables than the related CFLP. To 


apply cross decomposition to CFLPSS we first reformulate the problem as follows: 


ssi) W: i 


(ss2)  X;: j Vi 
Xj ij E{0,1} Vij 
52 >4, 4 

(SS 3) be i j 
y, €{0,1} Vi 


This formulation partitions the constraints into three sets: those that include only the 
facility location decisions (Y), those that include only customer assignment decisions (X,), and 
those that involve both types of decisions (W). Cross decomposition exploits this structure. 
Note that the constraint in CFLP that limits the assignment of customer demand to the 
available supply of open facilities (CFLP 2) has been replaced in CFLPSS by a constraint in 
X, that limits customer demand to the available supply of the facility. Thus the variable upper 


bounds (now in W) are no longer redundant. 
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A. GCD APPLIED TO CFLPSS (CX) 

Holmberg [1994] describes the problems and the general procedures to apply cross 
decomposition to CFLPSS. We adopt his notation with only minor modifications to label and 
classify the various decomposition problems. In this section we formulate the problems 
associated with generalized cross decomposition (GCD) and describe our implementation of 
GCD. We highlight where appropriate the implementation details not provided by Holmberg 
[1990, 1994]. In the next section we formulate the problems associated with one of the six 
alternate ordinary cross decomposition algorithms (CD6) described but not implemented by 
Holmberg and describe our implementation of that decomposition. 

1. Lagrangean Relaxation 

Our implementation of Lagrangean Relaxation (DW3) iterates between the 
subproblem and master problem described below. (Holmberg [1994] describes the 
subproblems and master problems for three possible dual decomposition algorithms that he 
labels Danzig-Wolfe decomposition.) When applied to CFLPSS this algorithm converges to 
the Lagrangean relaxation of the constraints in the set W. 

a. The Dual Subproblem (SSDS) 


(SSDS) 


maximize Dg (Av) +8 (V)- D4, 
i J 


YX 
Where for fixed dual values we solve the following subproblems: 
minimize aoe, ots 
2 (= 2 f; A yj 


y 3 j 
issps1) DSY: 


> 51°); > od, 
i J 


y, {0,1} Wi 


subject to: 
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minimize bor +A, 4+V; 5)°%; 


xeX, 7 


2A, 1) = 


issps2) DSX.: Vi 
t d : < 
pM Pi 
j 


Xi j E{0,} Vj 


subject to: 


The dual subproblem SSDS is the Lagrangean relaxation of the 
CFLPSS relative to the constraints in W. This relaxation creates a set of independent binary 
knapsack problems. DSY selects facilities to open. DSX; assigns customers to facility 7 
without regard to either the open facilities selected by DSY or the assignments made to other 
facilities. This relaxed problem provides a lower bound and a set of cuts for the master 
problem. | | 


b. The Dual Master Problem (SSDM) 


maximize baa a ae > 2 


(SSDM) 
v,A : 


subject to: 


(t,) (t,) (t,) 
(SSDM 1) 5) s (3 Cig +%j “A, +X; 5 Vij) Vit, 


j 
(y) (ty) (,) 
(SSDM2) O° S DF Yi es yj ; Vi Vty 
i 


i,j 


(SSDM3) Vj >Q Vi,J 


The dual master problem SSDM is obtained by adding dual cuts (¢, and 

t,) each time the respective dual subproblem is solved. This problem provides an upper 
bound and dual values for the subproblem. To ensure that SSDM 1s initially bounded, we 
include a set of initial cuts corresponding to all facilities open and all customers assigned to 


each facility. Holmberg [1990, 1994] does not suggest this. 
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2. Generalized Benders Decomposition 


Generalized Benders decomposition (GBD) iterates between the master 
problem and subproblem described below. When applied to CFLPSS, it converges to the 
convexification with respect to the assignment variables (CX). (The convexification with 
respect to a set of integer variables is the optimal solution over the feasible region defined by 


the convex hull of the integer feasible extreme points in the set.) This is the best lower bound 


that we can obtain using decomposition [Holmberg 1994]. 


a. The Primal Subproblem (PSL) 


(PSL) 


minimize BCA, v- yA, +f ‘y, 
j j 
x 


subject to: 


est) ¥,,20 Vij 


Where for fixed dual values we solve: 


(x);4 4 
I (A,Vv)= 
7 xeX, 7 


@sL2) PSDS;: y 
QoaKen Ss: 
j viprri 
subject to: 4 J 


x;; €{0,1} VJ 


The primal subproblem, PSL, is the Lagrangean relaxation of CFLPSS 
with respect to the constraints W after the facility locations have been fixed. Since this 
problem is a nonlinear mixed integer program, we do not solve it directly. It may be solved 
approximately using subgradient optimization, or optimally using Danzig-Wolfe 


decomposition. The master and subproblems used to solve PSL optimally are given below 


followed by brief description of the solution method. 
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(1)  Danzig-Wolfe Master Problem to solve PSL (DML) 
maximize > 53” ee 4, 
(DML) j | 

v,A 


subject to: 
5® < ; a ape 
(DML 1) SD) SY; YiVij 
i i,j 


k k k ; 
(ML 2) 6n < a, Ci +Xhy A; +x} , Vij) Vi,k 
j 


(@ML3) V; 2 0 Vi,j 
(2) Danzig-Wolfe Subproblem to Solve PSL (PSDS;) 


minimize Gt FAV; 5); 


(x) 79 4 
i (A,Vv) = : 
xeX, / 


PSDS;: Vi 
j 


x, €{0,1} Vj 


subject to: 


b. Solving PSL 
To solve PSL, we iterate between the Danzig-Wolfe subproblem, 
PSDS,, and the Danzig-Wolfe master problem, DML. PSDS, is identical to DSX;, from the 
dual subproblem (SSDM). DML is SSDM with the y cuts (SSDM 2) aggregated into a single 
constraint (DML 1). Holmberg [1994] provides these problems, but not the implementation 


details that follow in this section. These problems converge when 


gO AD2 5? Vi 
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To synthesize dual information for closed facilities, we use the following rule. 


; . 0 if ¢,,+4,20 

If y, =0 then v; , = ~ me 
lc, jt j] otherwise 

This rule ensures that (DML 2) is satisfied for the closed facilities. When DML is unbounded 

the facility choices are infeasible and we add the following constraint to PMCX and SSDS: 


y ¥21 where J = {i| y, =1}. 
ie 


Holmberg assumes feasibility for CFLPSS when the constraints are partitioned into the set 
X and the set Y. If this assumption does not hold the above constraint is required to achieve 
feasibility. We retain the cuts in DML that were found each time PSL was solved since these 
are valid at every iteration between PSL and PMCX (below). 

C. The Primal Master Problem (PMCX) 


minimize 
Py 


subject to: 


(t) t) t (t) (t) 
(PMCX 1) p= > (6, +4; +i) xt yA +>, ag Via; Vt 
i,j i j 


desi 2 Le) 


(@McCxX2) y, EY: i 
y; €{0,1} Vi 


The primal master problem PMCX is the convexification with respect 
to the assignment variables (x, ,). It is formed by adding a primal cut each time the primal 
subproblem is solved. It contains constraint (PMCX 2) to ensure that enough facilities are 
open to meet total demand. This is a necessary condition for primal feasibility. This 1s not 


a sufficient condition for primal feasibility; therefore, additional constraints are required. We 
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detect infeasibility when DML is unbounded in PSL and add a constraint (described in section 
b. above) to achieve feasibility. 
3. Generalized Cross Decomposition 

Generalized cross decomposition (GCD) iterates between PSL and SSDS until 
a convergence test fails or these problems converge in objective value. GCD converges to 
the convexification with respect to the assignment variables (CX). This is the better of the 
two lower bounds provided by DW3 and GBD if they differ. If the GBD lower bound 
exceeds the DW3 upper bound, we switch to GBD to find the better lower bound. The 
convergence tests for GCD are described below. 

a. The Primal Convergence Test (CTPCX) 

CTPCX uses the following cuts: 


Ce + AD + VO hx _ pa + > i= 2 ae -y, < Upper Bound 
i,j j i j 

If all of these cuts are satisfied then GCD continues by solving PSL. CTPCX uses the cuts 
in PMCX to determine if PSL can improve the upper bound. If any cut is not satisfied we 
solve PMCX. CTPCX also fails if PMCX has not been solved for a fixed number of 
iterations. 

b. The Dual Convergence Test (SSCTD) 
SSCTD uses the following cuts: 


Ci °%i + mh +d xi) — y,” a x,t) - 1$-2. > Lower Bound 
2 ia 4 j | 


ij 
If all cuts are satisfied then GCD continues by solving SSDS. SSCTD uses the cuts in SSDM 
to determine if SSDS can improve the lower bound. If any cut is not satisfied we solve © 
SSDM. 
B. CROSS DECOMPOSITION APPLIED TO CFLPSS (LX) 
We can relax the integer restrictions on the customer assignment variables and obtain 


another decomposition that converges to the same objective value as the linear programming 
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relaxation of CFLPSS with respect to the assignment variables. This form of decomposition 
is similar to the decomposition described by Van Roy [1986]. In this form we replace the 
primal master problem PMCX with PMLX and the primal subproblem PSL with PSLX. We 
~ also replace the primal convergence test CTPCX with CTPLX. The dual problems and the 
dual convergence test remain unchanged. These new primal problems and primal 
convergence test are detailed below. 
1. Benders Decomposition on LX 

Benders decomposition on LX (BD6) iterates between the following 
subproblem and master problem and converges to the linear programming relaxation of the 
assignment variables (LX). 

a. The Primal Subproblem (PSLX) 


- minimize > fi Ji +> 6, Xi ij 
ij 


(PSLX) 
subject to: 
yao) a 
@stxt) W i 
x; 4S); [v;,] Vis 
S45 xij $5; [a;] 
(PSLX2) X;: 4 J Vi 
x; ; 20 Vj 


The primal subproblem (PSLX) is the linear programming relaxation 
of the original problem (CFLPSS) with respect to the assignment variables after the facility 
locations have been fixed. It provides an upper bound and a set of dual values to form the 


cut for the master problem (below). 
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b. The Primal Master Problem (PMLX) 


minimize p 
» 
subject to: 


(PMLX1) 2 yay - 50), +), f-oe -y, Vt 
] i 


j i j 


dus W224; 
@MLX2) IY: i j 
y, €{0,]} Vi 


The primal master problem is obtained by adding a cut (PMLX 1) each 
time the PSLX is solved. It fixes the facility location variables for the next decomposition 
problem, and provides a lower bound. 

2; Cross Decomposition on LX 
Cross decomposition on LX (CD6) iterates between PSLX and SSDS in the 
subproblem phase until a convergence test fails or these problems converge in objective value. 
The new convergence test for GCD and procedures for objective convergence are described 
below. 
a. The Primal Convergence Test (CTPLX) 
CTPLX uses the following cuts: 


A _ Sy eo? +), fi- by vi -y; < Upper Bound 
j i i j " 

If these cuts are all satisfied then CD6 continues by solving PSLX. CTPLX uses the cuts in 

PMLX to determine if the upper bound can be improved. PMLX is solved if any cut 1s not 


satisfied or if PMLX has not been solved for a fixed number of iterations. 
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b. Convergence to a Bound | 
CD6 can converge to the better of the bounds provided by BD6 and 
DW3, even if it is not known which is better in advance. If during cross decomposition the 
current BD6 lower bound exceeds the current DW3 upper bound, then cross decomposition 
switches to BD6. Likewise, if the current DW3 lower bound exceeds the current BD6 upper 
bound, then the algorithm switches to DW3. 
C. DECOMPOSITION APPLIED TO A STATIONING MODEL 
In this section we describe the application of a decomposition algorithm to solve a 


simplified unit stationing model. That model is given below. 


(sm) minimize YS yt pe j *ij 
i i,j 
subject to: 


Yxy=l (al vy 
i 


qusM1) W: 
x,,<y; [vj] Vis 
Dik x, 55, ([@j,~] VE ER 
j 
(USM 2) X;: dik Xi ij SS; % [773.5%] Vik ER, Vi 
Sia Yi > dit [O,] Wk ER, 
qusM3) Y: 4 # j 


y, €{0,]} Vi 


This unit stationing model (USM) is an extension of the CFLPSS model that includes 
multiple resources. These resources are separated into two categories. One category (R,) 
is resources that are consumed by the presence of the unit such as motor pools and billets. 


The second category (R,) is resources such as ranges and maneuver acres that are required 
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by the unit, but are shared for the most part with other units. This simplified model captures 
the essential aspects of OSUB [Dell e#. al. 1994]. Note that this model has the same decision 
variables as the CFLPSS model with the same number of customers and facilities. 

We do not solve the model directly, instead we satisfy the constraints associated with 
the second resource category (R,) by variable reduction. That is, we fix the assignment 
variable to zero for any unit and facility pair that is not feasible due to a shortage of one of 
these resources. 

1. Benders Decomposition 

Benders decomposition on LX (BDLX) for the unit stationing model iterates 
between the following subproblem and master problem and converges to the linear 
programming relaxation of the assignment variables (LX). This decomposition is similar to 
BD6. 

a. The Primal Subproblem (USMPS) 


minimize 


(USM PSLX) DFID i iy 
Xx ; ra 
i,j 
subject to: 


Yaj=l (aj) VW 
i 


(USM PSLX1) W: 


x5), ([¥,j,] Vid 
Vidi Xj Six [ig] Wk ER, 


(USMPSLX2) X;1 4 OSX;,S1 Vj Vi 


x, ; =O if ok ER |(s,,/d; 4) <1 


The primal subproblem is the linear programming relaxation of USM 
after the facility location variables have been fixed. Note that the relaxed assignment variable 


is zero if the facility can not provide enough of a resource in the second resource category to 
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the unit. USMPS provides an upper bound on LX and a set of dual variables that form the 


cut for the master problem (below). 
b. The Primal Master Problem (USMPM) 


minimize 
(USM PMLX) 
Py 
subject to: 
(t) (f) (t 
i j i,j i,kER, 
(USM PMLX2) Y: i j 


y, €{0,1} Vi 


The primal master problem is formed by adding a cut (USM PMLX 
1) each time the primal subproblem is solved. Note that the constraints in USM PMLX 2 
require that the open facilities provide enough of each resource in the first category to meet 
the total demand. This is a necessary, but not a sufficient condition for feasibility. USMPM 


provides a lower bound and a set of facility locations for the primal subproblem. 
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IV. COMPUTATIONAL EXPERIENCE 


A. DESCRIPTION OF THE PROBLEMS 

Beasley [1990] provides a library of optimization problems that include capacitated 
facility location problems such as those used by Van Roy [1984] to test cross decomposition; 
however, many of these problems are infeasible for the CFLP with sole sourcing and the 
remaining problems are directly solved rapidly by CPLEX [1994]. We therefore generate 
more difficult random test problems of various sizes. 

Our problem generation scheme is to specify a range of values for the data and to then 
randomly select the data value from a uniform distribution over this range. We first specify 
the number of customers and facilities (the problem name follows the convention, e.g., 10C20 
is ten facilities and twenty customers) and the minimum and maximum customer demand. 
The minimum supply for each facility is equal to the average demand times the average 
number of customers per facility. The maximum supply is equal to three times the minimum 
supply. Arbitrarily, the minimum fixed cost for each facility is equal to the supply, the 
maximum fixed cost is equal to twice the supply, the minimum variable cost for each customer 
and facility is equal to the customer demand, and the maximum variable cost is equal to twice 
the customer demand. 

B. RESULTS AND ANALYSIS 

Tables 1 to 14 below present results from applying the decomposition algorithms 
described in the previous chapter to our test problems. Each table contains a measure of the 
problem size (number of facilities and customers), the results for each algorithm and some 
general comments about the problem. The results for each algorithm include the time needed 
in processor (CPU) seconds on an IBM System 6000 Model 590H, the lower bound 
determined, the number of times that the various problems (SSDM in the DM column, SSDS 
in the DS column, PMLX and PMCX in the PM column, and PSLX and PSL in the PS 
column) are solved and remarks about the algorithm's performance. We attempt to run each 
algorithm until it converges to within 0.01 (the upper bound minus the lower bound all 


divided by the lower bound) of the best lower bound obtained by the algorithm. We truncate 
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any algorithm that fails to converge after 7,200 CPU seconds or 250 iterations of any problem 
or 1,000 iterations of PSL. 

We use CPLEX [1994] to obtain the optimal objective function value of the linear 
programming (LP) relaxation of the problem. The lower bounds obtained by the 
decomposition algorithm may be slightly less than the lower bound from the LP relaxation due 
to the convergence and truncation criteria for the decomposition algorithms described above. 
After a sufficient number of iterations, the decomposition algorithms will determine bounds 
that are equal to or better than the objective function value of the LP relaxation. 
| We also attempt to solve the problems with the CPLEX [1994] mixed integer solver 
to provide a basis for comparison. This commercial solver employs state of the art branch 
and bound techniques that include preprocessing, heuristic rounding to obtain the first integer 
solution, cut-off and shortcut techniques, clique and cover cuts, and numerous other features. 
We report the objective function value of the optimal solution or the best feasible solution for 
each problem. When the optimal solution is not known, we report the quality of the branch 


and bound solution in the remarks. 
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owe [ en 1 _ 
GCD 272.91 2,643.0 5 13 (405 iterations in PSL) 
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Table 1. BD6 converges most rapidly and determines a good lower bound. CD6 performs well, but does 





not determine the best possible bound. GBD determines the best bound. GBD and GCD converge in a 
reasonable number of iterations, but their time performance is poor. PSL determines an upper bound in GBD 


and GCD that is equal to the optimal objective function value. All five decomposition algorithms find bounds 


that are better than the LP relaxation. 2 


PMLX & PSLX converge 


15 (839 iterations in PSL) 
PMCX & PSL converge 





Table 2. BD6 converges most rapidly and determines a good lower bound. CD6 performs well and 
determines the best lower bound. GBD and GCD converge in a reasonable number of iterations, but their time 
performance is poor. PSL determines an upper bound in GBD and GCD that is equal to the optimal objective 
value. All five decomposition algorithms find bounds that are better than the LP relaxation. CPLEX branch 


and bound performs well. 











Problem: 15C30 Facilities: 15 ae" 30 
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GCD 5,263.52 3,614.0 
4,613.95 3,641.0 


Table 3. BD6 determines the best lower bound and is the only algorithm to improve the bound from the 






(989 iterations in PSL) 
0.076 convergence gap 
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linear programming relaxation. CD6 does not obtain the same bound because PSLX and SSDS converge to 
within 0.01 first. DW3 converges to a bound in the least amount of time. PSL determines an upper bound in 
GBD and GCD that is equal to the optimal objective function value. GBD and GCD terminate before the 
problems converge to within 0.01 of the bound. 
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128.37 6,416.6 SSDM & PMLX converge 
CPLEX 6,091.50 6,960.0 Pe solution within 0.072 of optimal 


Table 4. BD6 converges rapidly to the best lower bound among the decomposition algorithms. None of the 
decomposition algorithms determine a bound that is better than the LP relaxation. This indicates that a 
convergence criteria of 0.01 truncates the algorithms too soon for this problem. GBD and GCD fail to converge 
or even provide reasonable bounds. CPLEX branch and bound obtains a feasible solution that is within 0.072 


of optimal. 








Problem: 20C40 Facilities: 20 Customers: 40 


Iterations 


LP Relaxation : 


SSDM & SSDS converge 


(635 iterations in PSL) 


0.072 convergence gap 





Table 5. BD6 determines the best lower bound. CD6 and DW3 converge more rapidly to a weaker bound. 


No algorithm obtains a bound that is better than the LP relaxation of the problem. GBD and GCD terminate 
before they converge to within 0.01 of the bound. 
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Table 6. BD6 converges rapidly to the best lower bound among the decomposition algorithms that converge. 


DW3 converges to the worse bound and requires more time that BD6 or CD6. CD6 converges to the same 
bound as DW3. If CD6 switches to BD6 when SSDS and SSDM converge instead of halting, it can detect the 
difference in the bounds and converge to the better of the two. 
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GBD 7284.57 5,406.0 (454 iterations in PSL) 





0.063 convergence gap 


Table 7. GCD determines the best lower bound. BD6, CD6 and DW3 converge more rapidly to a similar 


bound. No algorithm obtains a bound that is better than the LP relaxation of the problem. GBD halts before they 
converge to within 0.01 of bound. Neither GBD nor GCD determines a good lower bound. 
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Table 8. GBD determines the best lower bound. GCD converges to a similar bound in much less time. Every 


algorithm obtains a bound that is better than the LP relaxation of the problem. CD6 determines after two iterations 


which decomposition algorithm converges to the better bound and switches to DW3. 
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SSDS & PSLX converge 
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0.019 convergence gap 
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0.017 convergence gap 


CPLEX | 3,916.95 12,312.0 Be solution within 0.126 of optimal 


Table 9. GBD and GCD determine the best lower bounds, but do not converge. DW3 converges to a better 


bound in less time. BD6 and CD6 converge very rapidly to similar bounds. 
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0.10 convergence gap 


CPLEX 7,396.27 8,253.0 solution within 0.115 of optimal 


Table 10. GBD determines the best lower bound, but does not converge. GCD converges to a different bound 
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in much less time. DW3 converges to the same bound as GCD in even less time. BD6 and CD6 are extremely 


fast, but converge to weaker bounds. 
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Table 11. GBD determines the best lower bound, but does not converge. GCD converges to the same bound 


in slightly less time. DW3 converges to a good bound in much less time, but CD6 obtains the same bound even 
faster. CD6 again detects the best algorithm and switches to DW3 after two iterations. BD6 is extremely fast, 


but converges to the weakest bound. 


Problem: 30C75 (Facilities: 300 
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13,539.7 
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Table 12. DW3, CD6 and BD6 all converge to similar bounds. BD6 converges rapidly. CD6 performs well. 
DW3 does not perform as well. GBD and GCD fail to obtain usable bounds. 
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Table 13. DW3 determines the best lower bound. BD6 and CD6 converge more rapidly converged to a similar 
bound. GCD and GBD do not converge. 


Problem: 40C75 Facilities: 40 Customers: 75 
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(145 iterations in PSL) 
0.063 convergence gap 


0.082 convergence gap 
CPLEX 7,200.34 16,119.0 Pf solution within 0.114 of optimal 


Table 14. GBD and GCD determine the best lower bound, but do not converge. DW3 converges to a good 







bound, but is also very slow. CD6 and BD6 both obtain reasonable bounds very rapidly. 
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BD6 converges to a lower bound rapidly for all problems. DW3 is too slow to 
converge to a lower bound for some problems. DW3 performs poorly for larger problems 
and problems with smaller convergence tolerances because later iterations between SSDS and 
SSDM required significantly more time to solve than previous iterations. 

CD6 converges to a lower bound quickly for most problems. It converges rapidly 
enough that it is possible to alter the halting criteria to obtain a better bound. For example, 
CD6 can continue for some number of iterations beyond the convergence of PSLX and 
PMLX and possibly determine whether DW3 will converge to a better bound. 

GBD and GCD usually converge to a lower bound too slowly, but they often do 
produce good bounds in a reasonable number of iterations. GBD and GCD solve PSL by 
Danzig-Wolfe decomposition at every iteration and this requires too much time. 

C; USING THE LOWER BOUND 

Tables 15 to 18 below present results from using the decomposition problems to 
produce a feasible solution to the CFLPSS for some of the previous problems. To obtain this 
solution we use the CPLEX branch and bound algorithm to solve the problem obtained by 
fixing the facility decisions to the configuration last determined by the decomposition 
algorithm and setting the lower cut-off to the lower bound determined by the decomposition 
algorithm. For each algorithm the table contains the added time in seconds to obtain the 
solution, the total time in seconds, the objective function value, and the quality of the solution 


measured against the lower bound determined by that decomposition algorithm. 
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Problem: 10C20 





273.01 





Table 15. All algorithms determine feasible solutions with objective function values that 


are within at least 0.085 of the optimal objective function value in significantly less time 


than branch and bound. 
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2,590.76 


Table 16. All algorithms determine feasible solutions. The first problem solved by BD6 


is not feasible. CD6 determines a solutions with an objective function value that are within 


at least 0.03 of the optimal objective function value in less time than branch and bound. 
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Problem: 15C30 Facilities: 15 
Algorithm 








Remarks 






Table 17. GBD and GCD did not determine feasible solutions. All other algorithms did 


determine feasible solutions faster than branch and bound, but the quality of these solutions 
is not always good. 


6th problem 





Table 18. All three algorithms determine feasible solutions that were better than the best 


known branch and bound solution. CD6 performed significantly better than all other 
algorithms finding a solution with an optimal objective function value within 0.027 of 


optimal in less than three minutes. 


This procedure to find a feasible solution for the CFLPSS does not always work well 
for the larger problems in this study. Computational testing indicates that it is often necessary 
to solve the problem several times with different facility configurations to find a feasible 
solution and in many cases no feasible solution is found. A better approach may be to check 
for feasibility each time customers are assigned to facilities during the decomposition 


algorithm and retain feasible facility configurations. 
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D. A GRAPHICAL EXAMPLE OF CONVERGENCE 

In this section we graphically compare the convergence of cross decomposition (CD6) 
with Benders decomposition (BD6) and Danzig-Wolfe decomposition (DW3) for problem 
20C40 shown in Table 5. For the problems tested the convergence characteristics of these 
decomposition methods on this problem are typical. In figures two through four below, the 
y-axis is labeled with objective function value and the x-axis is labeled with the CPU time at 
the completion of each iteration of the decomposition algorithm. We plot the best bound 


determined by the decomposition problem at each iteration. 


Bound Value 





0.21 0.31 0.43 0.54 0.63 0.75 0.88 1 
Time (cpu seconds) 


Figure 2. We solve PMLX and PSLX at every iteration. PSLX rapidly determines a tight upper bound 
and does not improve this bound after the second iteration. PMLX gathers cuts and improves the lower 
bound steadily at every iteration. BD6 halts when the objective function value of PMLX is within 0.01 of the 
upper bound determined by PSLX. 
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0.2 1.38 3.79 10.18 31.07 46.84 53.62 
Time (cpu seconds) 


Figure 3. DW3 solves SSDM and SSDS at every iteration. These problems do not provide good bounds 


initially, but after thirty-five iterations without progress they begin to converge. Note that the time to solve 
each iteration increases markedly as the decomposition algorithm progresses. The first thirty iterations 
require less than half as much time as the next ten iterations. The last sixteen iterations account for over forty 
percent of the total solution time. SSDS determines a good upper bound after forty iterations and only 
improves slightly after that point. SSDM gathers cuts and improves steadily after forty iterations. DW3 
halts when the objective function value of SSDM is within 0.01 of the best upper bound determined by 
SSDS. 
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0.11 0.31 0.52 0.7 1.34149 1.7 1.84 2.16 234 2.74 2.99 3.31 3.52 3.81 4.03 4.36 4.73 5.07 5.74 6.08 648 6.81 
Time (cpu seconds) 


Figure 4. CD6 solves SSDS and PSLX at every iteration. These decomposition problems determine 
moderately good bounds initially, but do not converge rapidly. PSLX never improves the upper bound after 
the second iteration. When a convergence test fails, a master problem is solved. SSDM is solved more 
frequently than PMLX, but this problem never provides a good upper bound. PMLX provides a good lower 


bound when first solved on the 12th iteration. CD6 halts when the objective function value of PMLX 
converges with the upper bound determined by PSLX. 


It is possible that the lower bound determined by DW3 is superior to the bound 
determined by CD6 for this problem. To determine if this is the case we would iterate 


between SSDM and SSDS after PSLX and PMLX converged until either they converged or 
the bound from PMLX exceeded the bound from SSDM. 
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E. | COMPUATATIONAL EXPERIENCE WITH USM 
1. Description of Problems 

We generate random test problems of various sizes for the unit stationing 
problem. The problem generation scheme is similar to that described previously for CFLPSS 
with slight modifications for multiple resources. We first specify the number of units and 
facilities (the naming convention is e.g., 24U36 is 24 facilities and 36 units) and the minimum 
and maximum unit demand. The minimum supply of resources of type R, (those resources 
not shared by units) for each facility is equal to the average demand times the average number 
of units per facility. The maximum supply of resources of type R, is equal to three times the 
minimum supply. The minimum supply of resources of type R, (those resources shared by 
units) for each facility is equal to the minimum demand. The maximum supply of resources 
of type R, is equal to the maximum demand. Arbitrarily, the minimum fixed cost for each 
facility is equal to the maximum supply among the resource categories, the maximum fixed 
cost is equal to twice the minimum fixed cost, the minimum variable cost for each unit and 
facility is equal to the unit demand, and the maximum variable cost is equal to twice the unit 
demand. 

2. Results 

Table 19 below presents results from the Benders decomposition algorithm 
(BD6) applied to the unit stationing model, the branch and bound algorithm applied to the 
problem obtained by fixing the facility configuration to that of the last iteration of BD6, the 
branch and bound algorithm applied to the unit stationing model (USM), and the simplex 
algorithm applied to the linear programming relaxation of the unit stationing model. For each 
problem the table contains for BD6 the total time in seconds, the lower bound, and the 
number of iterations for BD6 to converge; for the problem to assign units to the facilities 
selected by BD6 the total time in seconds and the objective function value; for branch and 
bound the total time in seconds and the objective function value; and for the linear 


programming relaxation the objective function value. 
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Problem . BD6 


— indicates that no feasible integer solution is found. 





Table 19. BD6 converges rapidly and determines a good lower bound for all problems 


except 18U36. BD6 finds bounds for all problems that are better than the LP relaxation. BD6 
finds feasible integer solutions for all problems except 24U36. All CPLEX branch and bound 


solutions are the best known integer solution. 


These results indicate that a Benders decomposition may be a good alternative to 


branch and bound for finding good feasible solutions to USM. 
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V. CONCLUSIONS 


A. PERFORMANCE OF THE DECOMPOSITION ALGORITHMS 

Benders decomposition (BD6) is the best algorithm, especially for larger problems, 
because it converges to a lower bound most rapidly. Lagrangean relaxation (DW3) performs 
well for some small problems, but later iterations of this algorithm require significantly more 
time to solve than the earlier iterations. Cross decomposition (CD6) is a good algorithm for 
medium sized problems. It converges rapidly enough that it may be possible to continue the 
algorithm until it is clear whether BD6 or DW3 will converge to the better bound. 

Generalized Benders decomposition (GBD) and generalized cross decomposition 
(GCD) are too slow, but they should not be abandoned completely. They are guaranteed to 
converge to the best possible lower bound and they do converge in a reasonable number of 
iterations. The primary problem with GCD and GBD is the poor performance of PSL. The 
next section discusses techniques to address this shortcoming. 
B. TOPICS FOR FURTHER RESEARCH 

1. Determining Upper Bounds 

The decomposition algorithms presented in this thesis find lower bounds on 

the optimal objective function value. In the process they may also identify a feasible set of 
facilities. The lower bound is useful to judge the quality of a feasible solution. Fixing the 
facility variables determined in the last iteration of the decomposition algorithm and solving 
the resulting problem is typically an efficient method to obtain a feasible solution. This 
method fails if the set of facilities does not contain a feasible set of customer assignments. 
Other methods for obtain upper bounds and feasible solutions should be explored. 

2. Improving the Efficiency of the Algorithms 

It is possible to improve the efficiency of the current implementations in 

several ways. For the decomposition algorithms considered in this thesis there are techniques 
to strengthen the cuts currently determined in the subproblems. Application of these 
techniques would reduce the number of iterations required to achieve convergence and may 


prove computationaly attractive. 
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As implemented for this thesis, GBD and GCD are not viable solution 
methods. However, they are theoretically very attractive because they converge to the best 
bound that can be obtained from the decomposition algorithms we considered. In GBD and 
GCD the most expensive problem is PSL. Three alternate solution approaches for PSL 
should be considered. 

(1)  Useaheuristic method (such as solving a relaxed problem and rounding to an 
integer solution) to solve PSDS, until PSDS; and DML near convergence, then 
switch to an optimal solution method. 

(2) Solve the binary knapsack problems in PSL with a pseudo-polynomial 
algorithm (such as dynamic programming). These methods are often superior 
to branch and bound. 

(3) Solve PSL using subgradient optimization instead of Danzig-Wolfe 
decomposition. 

2. Implementing the Stationing Model 
The motivation for this thesis research was the unit stationing model. The 

Benders decomposition method for the stationing model outlined in this thesis can be 
extended to cross decomposition by implementing a Lagrangean relaxation algorithm. 
Computational experiments should be performed with these decomposition algorithms on a 
larger set of problems. 

4, Improving the Experimental Design 
The computational results reported in this thesis are for a limited number of 

instances of each problem size and for one set of parameters that generate the problems. It 
would be better to test the algorithms against a random sample of problem instances for each 
problem size. It would also be useful to generate problems with different sets of parameters. 
This may lead to additional work to establish a set of metrics for CFLP problems and to 


quantify the difficulty of problems in terms of these metrics. 
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APPENDIX. BENDERS DECOMPOSITION 


In this section we derive the Benders master problem for the capacitated facility 
location problem (CFLP). In Benders decomposition we identify complicating variables that 
can be fixed to produce a simpler problem. For the CFLP the binary variables (y) associated 
with the facility location decisions complicate the problem. With the y variables fixed the 
CFLP is a simple linear program: a bipartite network with gains that provides a minimum 
cost assignment of customer demand to facilities. 

The problem below is a reformulation of the CFLP with an outer optimization over 


the complicating variables (y), and an inner optimization over the simple variables (x). 


minimize 
x 


os “Xi j 
i,j 


subject to: > 4%) =! [A,] W 
DAF od 
yy 1 S51: [4;] Vi 
j 
xij SY; [V; ] Vi, J 
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i j 


y,; €{0,1} Vi 


The added constraint in the outer optimization ensures that the facility configuration 1s 


feasible. 
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With the facility configuration fixed we replace the inner optimization with its linear 


programming dual: 


maximize Dy = >» Yi ki ->y, ey 
j td 
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minimize 
Df yit 

y j 
subjectto: A,-d;-H;-V;j S$, [x,;] Vis 


w,20 Vi v,,20 Vij 


subject to: Dee -y,2 yd j 
j 


y; €{0,1} Vi 
Let 


T= {t | (a®, a”, vy) is a (dual) feasible extreme point solution} 


be the index set of all (dual) feasible extreme point solutions of the inner optimization problem 


above. Then the previous problem can be written as 


minimize |maximize 

(t) (t) (t) 
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This problem is equivalent to the original problem since the dual of the inner optimization 
problem attains its optimal solution at one of a finite number of extreme points. (The dual 
problem is bounded if the primal problem is feasible.) The problem above can also be written 


as the mixed integer program below. 


minimize 
p 
Py 


subject to: p2 qs fi -S:° uw” a yk y, WteT 
J I i i iJ i 
j i j 


Desi? Dd, 
ij 


y, €{0,]} Vi 


The Benders master problem is a relaxation of this problem obtained when only a subset of 
the constraints associated with the index set J are known. The Benders master problem 
provides a feasible facility configuration and a lower bound on the optimal objective function 
value of the original problem. 

The set 7 is usually large and only a subset of these constraints are binding in the 
optimal solution. The strategy adopted by Benders decomposition is to solve a relaxed 
master problem that contains a subset of these constraints to obtain a feasible facility 
configuration and to then solve another problem, the Benders subproblem, to determine if a 
lower cost facility configuration exists. Ifthe facility configuration from the relaxed Benders 
master problem is not optimal, the subproblem provides a new constraint that is violated by 
this facility configuration. This new constraint, called a Benders cut, is one of the constraints 
in the index set 7 that is not already in the master problem. The problem below is the 


Benders subproblem. 


f 
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This problem is identical to the inner optimization problem in the reformulation of the CFLP. 
The Benders subproblem is a restriction of the CFLP obtained by fixing the facility locations. 
The optimal solution to the Benders subproblem provides a set of dual variable values that 
form a Benders cut and an upper bound on the optimal objective function value of the original 
problem. The objective function value of the Benders subproblem is equal to the optimal 
solution of the original problem when the facility configuration is optimal. 

Benders decomposition starts with a feasible facility configuration and iterates 
between the Benders subproblem and the Benders master problem. At each iteration, if the 
facility configuration is not optimal, the subproblem provides a new Benders cut. This cut 
is a violated constraint from the constraints associated with the index set 7. In the extreme 
case, after a finite number of iterations the subproblem produces all of the constraints 
associated with the index set 7 and the master problem is equivalent to the original problem. 
Thus, after a finite number of iterations the master problem and the subproblem must 
converge in objective function value to the optimal objective function value of the original 


problem (CFLP). 
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